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Abstract 

In this paper we utilize the covariance of Riccati equation with respect to linear fractional 
transformations to define classes of conformally equivalent second order differential equations. 
This motivates then the introduction of fractional Darboux transformations which can be 
recognized also as generalized Cole-Hopf transformations. We apply these transformations 
to find Schrodinger equations with isospectral potentials and to the linearization of some new 
classes of nonlinear partial differential equations. 



1 Introduction 



Darboux transformations [1,2,3] of the form 



if, 



A(x) + B{x) 



d_ 

dx 



(f){x) 



1.11 



have many different applications in mathematical physics [9,10,11,12,13] and the theory of dif- 
ferential equations. Among these one can find the Factorization method (raising and lowering 
operators for one or several coupled differential equations [1,4]) separation of coupled Schrodinger 
equations [6] and applications to nonlinear integrable equations. Furthermore with proper gen- 
eralization this transformation has been applied to multidimensional problems [6,10,14,16] and 
discrete systems [15,17]. 

In its simplest context the transformation ( |1 . 1| ) "connects" the solutions of two Schrodinger 
equations with different potentials u(x),v(x) i.e. 



0" = ( u (x) + A)0 



1.2) 



^" = ( v (x) + X)ip. 



1.3) 



Letting B(x) = 1 one can easily show that in order for eqs ( 1.2|) , (|1.3| ) to be related by the 
transformation ( |1.1[ ) A(x),u(x),v(x) must satisfy; 



A" + v! + A{u - v) = 



(1.4) 



2A' + u - v = 



Eliminating {u — v) between these equations and integration yields 



1.5) 



A' - A 2 + u = -v 



1.6) 



where v is an integration constant. Eq. (|1.6j ) is a Riccati equation which can be linearized by the 
transformation A = — C'/C which leads to 



C" = (u(x) + 

1 



(1.7) 



Thus ( is an eigenfunction of the original eq. ( |1.2| ) with A = v. From ( |1.5| ) we now infer that 

v = u-2{lnQ" (1.8) 
i.e. a Darboux transformation changes the potential function u(x) by Au = —2(lnQ" where ( is 



an arbitrary eigenfunction of (|1.2| ). 

Our objective in this paper is to introduce fractional Darboux transformations which are defined 
as transformations of the form 

A(x)4>(x) + B(xW 
W C{x)<f)(x) + D{x)<j)' 1 ' 

and elaborate on some of their applications to linear and nonlinear differential equations. To 
motivate the introduction of these transformations we use the covariance of the Riccati equation 
with respect to linear fractional transformations. We then show that this covariance induces 
an equivalence relation between different second order differential equations whose solutions are 
related by a transformation of the form ( |1.9| ) (Sec. 2). In Section 3 we consider the classification 
problem for these "conformally equivalent" equations. In Section 4 and 5 we explore then the 
application of the transformation ( |1.2| ) to Schrodinger equations and nonlinear equations. We end 
with some observations and conclusions in Section 6. 

2 Conformally Equivalent Equations 

Consider a linear second order differential equation 

p(x)w" + q(x)w' + r(x)w = (2.1) 

(with smooth enough coefficients). 
Introducing 

V= U - (2.2) 

w 

eq. ( P7T| ) is transformed into a Riccati equation 



However it is well known that Riccati equation is "covariant" with respect to the linear fractional 
transformation [5,7] 

= a(x)z(x)+i(x) 
V (3{x)z{x) + 8{x) 1 ' ' 

i.e. such a transformation takes a Riccati equation into another one in its class 

z' = F(x)z 2 + G(x)z + H{x) (2.5) 
In particular for eq. ( |2.3j ) we have; 



f =-v + (M) + 01 W>) + {a ' 13 - a0m (2 ' 6) 

G = -[2« 7 + (a5 + /3 7 ) f + 2/55 + (a'5 - a5') + ( 7 '/3 - 7/3')]/A (2.7) 



v p(a;) 

p(x) / ' Ip(a^) 



= _[ 7 2 + 7 5 f 441 + 5 2 f44l + (V* - lS')]/A (2.1 



where A = a5 — /3 7 7^ 0. 

However we can transform eq. ( |2.5| ) back to a second order linear differential equation by the 
transformation 

-n'Qr) 

to obtain 

u" - ( G(x) + u' + H(x)F(x)u = 0. (2.10) 

V F ( X )J 

We conclude then that the sequence of transformations ( |2.2|) , (|2.4j ), ( |2.9|) establishes an equivalence 
relation between linear second order differential equations and we shall say that such equations 
are "conformally equivalent" . 

To see how the solutions of eqs (|2.2j), ( ^.1(J| ) are related we note that from (P^2"|), ( |SO| ) we have 

Sw' — jw 



aw — f3w' 

and hence using ( |2T9| ) 

ii.(t.) = pxn I / P(t\ 

aw — (3w' 



(2.11) 



u(:r) = exp I / F(a;) — — ^dx] . (2.12) 
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We conclude then that up to exponentiation and integration the relationship between the solutions 
is given by a Fractional Darboux transformation. 

A(x)w + B(x)w' 
K ' C(x)w + D(x)w' 



We now explore some properties of this relationship when a,/3,j,S is eq. ( |2.4|) are constants. 
Under this restriction it is well known that the transformation (|2.4| ) is equivalent to a sequence of 
three transformations 

y = a 1 z 1 + b 1 (2.14a) 
z x = — (2.14b) 
Z2 = a 2 Z3 + b 2 (2.14c) 



Since affine transformations ( |2.14a| ), ( |2. 14cj ) induce only a minor change in the original equation 



(they just add constant to r(x)) it is interesting to find those equations which remain invariant 
with respect to the transformation y = -. To this end we let p(x), q(x), r(x) be polynomials in x 



P( x ) = ^2PkX k , q{x) = J2 VkX k , r(x) = ^r k x k (2.15) 

k=0 k=0 k=Q 



and solve the equations 

^ = H(x)F(x) (2.17) 

for the coefficients p k , qk, hk- When n = 2 we find that up to translations in x (i.e. x = x + a) 
there are only three equations with this property 

(1) 

w " - n 2 w = (2.18) 

(2) Chebychev equation 

(p - x 2 )w" - xw' + h w = (2.19) 
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(3) 



xx - 2 — w H w H x u> = (2.20) 

V Po/ Po Po 



Letting p — 1, g = § in ( |2.20| ) we have 

x(l — x)w/' — — h x 2 w = 

This equation does not appear in Kamke's book [5]. However it belongs to the class of Hill-Mathieu 
equations through the transformation x = cos 9. 

In Table 1 we present the effect of the transformation y = - on several classes of equations 
whose solutions are the Special Functions of Mathematical Physics. We would like to emphasize 
however that in general a, (3, 7, S in eq. (|2.4j) can be functions of x rather than just constants and 
consequently the conformal equivalence class of each of these equations can be very large indeed. 

Finally for the Schrodinger equation ( |1.2[ ) we find that when a, (3, 7, 5 are constant with a5 — 
/?7 = 1, Q(x), R(x) in eq. ( p.lO|) take the following form; 



3 2 u'(x) 

Q(x) = 2a 7 - 2f35(u(x) + A) - — ^ 1 >- - (2.21) 
R(x) = [5 2 (u + A) - 7 2 p 2 (w + A) - a 2 } (2.22) 



3 The Classification Problem 

To use the algorithm described in the previous section in a practical manner one must be able to 
determine when a differential equation 

w'l + qi(x)w[ + r\(x)wi = (3.1) 

is conformally equivalent to eq. (|2.1|). To this end one must solve the equations 

qi (x) = G(x) + (3.2) 

n(x) = H(x)F(x) (3.3) 
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for appropriate a(x), (3(x),j(x),8(x) subject to the condition 

To make progress towards the solution of this classification problem we shall consider two 
separate cases; /3(x) ^ and (3{x) = 0. 

A. p(x) ± 0. 

In this case we can assume without loss of generality that (3{x) = 1. The resulting expression 
for qi(x) can be used to solve for j(x) in terms of qi(x), a(x), 5(x) and the coefficients of eq. (2.1). 
Substituting this expression for j(x) in eq. (|4~3l ) (with /3(x) = 1) we obtain; 



n(x) = ^y- + --^ + R^x) (3.4) 

where Ri(x) which depends only on a(x) and the coefficients of eq. (|2.1|) is given by eq. (A.l). 
(To simplify this equation we let p(x) = 1). 

Eq. (|3.4j ) is a single ODE for a(x) and if a solution to this equation exists then equations 
( |2.1| ), ( p.l| ) are conformally equivalent. However in general eq. Q3.4| ) is nonlinear (and intractable) 
equation for a(x). Due to this circumstance it is appropriate to use eq. (3~4") in a "reverse" manner 



viz. assume a fixed functional form for a(x) and determine the corresponding form of Ri(x). This 
expression can be evaluated for various classes of 2nd order ODEs and an appropriate table can 
provide a quick reference to determine if a given equation is conformally equivalent to any of the 
classes that appear in the table (under the restrictions imposed on a(x)). 

Table 2 presents the different R\(x) that correspond to some classes of special functions when 
a(x) = 0. 

From this table we gather the following: 
Corollary 3.1 1. u" + uo = is conformally equivalent to Bessel equation of order 1/2 

2. Bessel equations of order 0, 1 are conformally equivalent to each other. 
B. p(x) = 

As can be expected this is a simple case to treat. 



6 



To begin with we can assume without loss of generality that S(x) = l(a(x) ^ 0) and obtain 
for q%(x), ri(x) the following expressions (we let p(x) = 1) 

qi(x) = 2a(x) + q(x) (3.5) 

o ()(y. yXJ 
r\{x) = a (x) + a(x)q(x) + r(x) H — . (3.6) 

Solving ( |3.5| ) for a(x) and substituting in ( |3.6| ) we obtain eq. ( |3.4| ) with Ri(x) given by 

Ri(x) = —~2~dx~ (X) ( ) 

(observe that i?(x) is independent of a(x)). 

To recast this algorithm in perspective we note that Forbenius method expresses the solution of 
(|2.1|) in terms of powers of x. The algorithms described above enable us to approach this problem 



differently by asking whether the solution of (|3.1| ) can be expressed naturally in terms of other 



classes of function (e.g. special functions) i.e. when (|3.1| ) is conformally equivalent to an equation 
whose solution is already known. 



4 Application to Schrodinger Equation 

In this section we consider the application of fractional Darboux transformations in the form 

A{x)<j>{x) + <}>'{x) 
W[X> B{x)(j){x) + 4>'{x) 1 ' 

to relate the solutions of equations ( |1.2|) , Differentiating ( |4.1| ) twice and substituting in 

( |1.3| ) using ( |1 . 2| ) to eliminate the higher order derivatives of <p(x) we obtain a polynomial which 

contains the monomials 3 , <p 2 <p', <f){4>') 2 and (</>') 3 . To annul this polynomial we let the coefficient 

du 

of each of these monomials be zero. By using the last three equations to express X,v(x) and — 

(JjUj 

and simplifying we obtain from the first equation (coefficient of 3 ) 

[(B - A)' - (B 2 - A 2 )]' - 2B[(B - A)' - (B 2 - A 2 )] = 0. (4.2) 
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This can be solved by the Ansatz 



[B - A)' - (B 2 - A 2 ) = (4.3) 



substituting this relation in the equations for A, v and — yield the following equations 

ax 

u(x) = A 2 -A' - c (4.4) 

v{x) = 2A(A -B)-c (4.5) 

A = c (4.6) 

where c is a constant of integration. Thus 

Au = v(x) -u(x) = A(A-2B) +A'. (4.7) 

To use these equations in practical applications we start with eq. (4.4) to find A. This Riccati 
equation can be converted to 

Ci = (u(x) + c)Ci (4.8) 



A=-$. (4.9) 
Ci 



by the transformation 



It follows then that is a solution of eq. (|1.2|) with A = c. Similarly if we substitute for A in 
fl4~3|) using (U) we obtain 

B' -B 2 = -u- c. (4.10) 

Thus (2 (where B = — — j is also a solution of ( |1.2[ ) with the same eigenvalue. It follows then 
that 

Au = [(/<i)f - 2(ln( 1 )'(ln( 2 )' - (ZnCi)" (4.11) 
which is completely different from ( |1.8| ) for Darboux transformation ( |1 . 1| ) . 
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5 Generalization of Cole-Hopf Transformation 

It is well known that Cole-Hopf transformation [8] 



(j)(x,t) 



(5.1; 



can be used to linearize the Burger's equation 
and transform it to the heat equation 

4>t = vcf) xx . 

It is now easy to recognize that the fractional Darboux transformation 

A(x)<j>(x,t) + C(x)(f> x (x,t) N 



*p(x,t) 



B(x)<f>(x,t) + D(x)(f> x (x,t) E 



(5.2) 



(5.3) 



(5.4) 



is actually a generalization of (5TT). It is natural therefore to apply it to the linearization of 
nonlinear partial differential equations. 

In the following we restrict ourselves to the classification of nonlinear equations which can be 
transformed to the heat equation through the transformation (|5.4|) . 

To this end we compute 

dip d 2 tp 
dt dx 2 

using ( |5.4|) and simplify the expressions thus obtained by replacing <f> t by (p xx at each step (we set 
v = 1 in eq. ( |5.3|) ). 

After some algebra we obtain 



A" + (2A' + C") ( 


4> J 


+ c 


( 4>xx \ 

\ <t> ). 


B + D 


if. 


) 



(5.5) 



B" + (2B' + D") ( 






B+D 


(t 


) 



E x 
~E 



However from 



Bip - A 
C-D-4)' 
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(5.6) 



Substituting^]^) and using the identity 

% = ( ^] + { Z£ ] . (5.7) 



. 2 

J x 



in eq. (|5.5| ) we obtain for the right hand side an expression which contains monomials of ip,ip x 
up to the fourth and second degree respectively. However if we restrict ourselves to second order 
nonlinearity it is enough to let D = and set B(x) = 1. The simplified form of eq. ( |5.5|) in this 
case is 

C\ip t - il> xx ) = -Cip 2 + [C'(2A + C) - C(C - 2A)' + 2Cijj x )iP 
- C(C + 2A)tfj x + AC(C + 2A)' - A{C" + C'A) 
+ C(C'A' - A"C). (5.8) 



In particular if A, B,C, D are constants eq. (|5.5|) simplifies and we obtain 

VH ~ Wxx = Zipx Bi) -C ' ' 

which is a generalization of Burger's equation. 

Similarly we can use the transformation (5.1) (with A,B,C,D being functions of x only) to 
relate a linear second order equation of the form ( [2.1|) (with w replaced by <p) to second order 
nonlinear ODEs. If we restrict ourselves to nonlinearities up to the third order we find that we 
must set D(x) = (and without loss generality let B(x) = 1). The general form of the resulting 
nonlinear second order equation is 

CV = 2ip 3 -(3qC + 6A + 2 C')^j 2 

+ [(-2 r + q' + q 2 )C 2 + (C" + 2 q C + 6 Aq)C 

+ QA 2 + 4C'A}^ + (r' + qr)C 3 (5.10) 

+ [(-?' - q 2 + 2 r)A + A" + 2 C'r]C 2 

+ [-3 gA 2 - (C" + 2 gC")^]C - 2 A 3 - 2 C'A 2 

Similar algorithm can be used to find other nonlinear equations that can be transformed to a 



given linear equation via a transformation of the form ( |5.4| ). 
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6 Conclusions 



We introduced in this paper Fractional Darboux transformations and elaborated on some of their 
applications. In view of the wide range of applications that regular Darboux transformations 
have in the literature we expect that this new transformation will generalize many known results. 
In particular the systematic adaption of this new transformation to the computation of Laddar 
operators for new potentials in the spirit of [1,4] has yet to be worked out. We hope to elaborate 
on these subjects in subsequent publications. 
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Appendix 



In this appendix we give the general form of Ri(x) in eq. (3.4). 



„ . . Ni(x) . . _ s 

«■(*) = ^ (A. 1) 



where 



JVi = 2(a 2 + aq + a' + r)a" - 3(a") 2 - 6(2a + (qa)' + r')a" 
+ 12(a') 3 + 6(4r + q' - q 2 )(a') 2 + [4(4r - 2q' - q 2 )a 2 
+ 2(g" + 8qr - 8/ - 2q 3 )a + 8r(2r + g') - Aq(qr + 2r') + 2r"\ol 
+ (4r - 2g' - g 2 )a 4 + 2(g" - 2/ - g 3 + 4gr - g'g)ct 3 (A. 2) 

+ [8r 2 - g 4 + 2g 2 r + [2q" - Qr')q - 3q" + 2r"]a 2 + [8gr 2 
+ 2(g" + g'g - 2r' - g 3 )r + 2r"g - 2gV - 6r'q']a 
+ 4r 3 + (2q' - q 2 )r 2 + (2r" - 2gr')r - 3(r') 2 

Di = 4[a(a + g) + a' + r] (A. 3) 

For brevity we suppressed the dependence on x in eqs. (A. 2), (A. 3). 
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Table 1: The effect of the sequence of transformations (2.2), (2.4), (2.9) on the equations in the 
second column when (2.4) is given by y = \ 

Class of equations Original equation Transformed equation 



Constant Coefficient 



Legendre 



Hermite 



Bessel 



Leguerre 



Chebyshev 



Hypergeometric 



au" + ku' + au = 

(1 - x 2 )w" - 2xw' + n(n + l)w = 
w" - 2xw' + 2nw = 
x 2 w" + xw' + (x 2 — n 2 )w = 
xw" + (1 — x)w' + nw = 
(1 - x 2 )uj" - xuj' + h 2 u = 
x(l - x)w" + [7 - (a + (3 + l)x]w' 
—a/3w = 



Same 

(1 - x 2 )u" + n(n + l)u = 
u" + 2xu' + 2nu = 
u" - f 2 +n \ v! + = 

u" + u' + ^u = 

X 

Same 

x{\ - x)u" + [l-7 + (oi! + /3 - l)x]% 
—a/3u = 
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Table 2: Form of Ri(x) in eq. (3.4) for some classes of special functions when a(x) = 







Constant coefficients 




c~\b 2 




Legendre 




n(n+ 1) 
1 -x 2 




Hermite 




-x 2 + 2n-l 




Bessel 


x 2 [4x A — 


3x 2 (4n 2 + 1) + 2n 2 (6n 2 - 5)] + 


n 4 (4n 2 + 1) 




4x 2 (x 2 — n 2 ) 2 




Leguerre 




n 1 

x 4 




Chebyshev 




(4n 2 - l)x 2 - 2(1 + 2n 2 ) 






A{x 2 - l) 2 n 4 




Hypergeometric 


[1 - (a-P) 2 ]x 2 


+ 2[(a + /? + l) 7 - (2a + l)P- 


(l + a)]x + 7 2 -l 


x 2 [x — l) 2 
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